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ABSTRACT 

The low-inclination component of the Classical Kuiper Belt is host to a population of extremely 
widely-separated binaries. These systems are similar to other Trans-Neptunian binaries (TNBs) in that 
the primary and secondary components of each system are of roughly equal size. We have performed 
an astrometric monitoring campaign of a sample of seven wide-separation, long-period TNBs and 
present the first-ever well-characterized mutual orbits for each system. The sample contains the 
most eccentric (2006 CHgg, e„i = 0.9) and the most widely-separated, weakly bound (2001 QW322, 
a/ Rh — 0.22) binary minor planets known, and also contains the system with lowest-measured mass 
of any TNB (2000 CF105, Mgys — 1.85 x 10^^ kg). Four systems orbit in a prograde sense, and three in 
a retrograde sense. They have a different mutual inclination distribution compared to all other TNBs, 
preferring low mutual-inclination orbits. These systems have geometric r-band albedos in the range 
of 0.09 — 0.3, consistent with radiometric albedo estimates for larger solitary low-inclination Classical 
Kuiper Belt objects, and we limit the plausible distribution of albedos in this region of the Kuiper Belt. 
Wc find that gravitational collapse binary formation models produce a similar orbital distribution to 
that currently observed, which along with a confluence of other factors supports formation of the 
cold Classical Kuiper Belt in situ through relatively rapid gravitational collapse rather than slow 
hierarchical accretion. We show that these binary systems are sensitive to disruption via collisions, 
and their existence suggests that the size distribution of TNOs at small sizes remains relatively shallow. 
Subject headings: Kuiper belt: general — planets and satellites: dynamical evolution and stability 



1. INTRODUCTION 

Most, if not all, minor planet populations are host to 
multiple systems, including binaries and trinaries (see 
Walsh 2009 for review). These multiples have a vast ar- 
ray of properties, from extremely short-period, contact 
binaries, to systems that are exceedingly widely sepa- 
rated and have mutual periods of many years. Some 
have tiny satellites, thought to be coUisional fragments 
blown off of their parent, and others have components 
of near-equal size and mass. The characteristics of these 
systems represent a treasure trove of information about 
the properties of the objects that compose them, the en- 
vironment they are embedded in, and the dynamical his- 
tory of their parent population. Through their orbital 
separations and periods, binaries offer the only way to 
measure the mass of these distant objects, which when 
combined with radius measurements determine these ob- 
jects' bulk densities — which in turn provide information 
about composition and physical structure (such as poros- 
ity). 

In the last 10 years, Trans-Neptunian populations have 
been found to host to a very high fraction of binary sys- 
tems. The binary fraction varies in sub-populations from 
~ 30% in the low-inclination component of the Classical 
Kuiper Belt to just a few percent in other dynamical 
classes (Noll et al. 2008a). Given the low interaction 
rates of the Kuiper Belt populations today, forming such 
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a large number of binary systems has proven a theoretical 
challenge, especially with the limited information avail- 
able for the components of these systems (eg., Goldrcich 
et al. 2002, Weidenschilling 2002, Noll et al. 2008b, 
Schlichting & Sari 2008a & b). 

Trans-Neptunian Binaries (TNBs) are distinguished 
from binary systems elsewhere in the solar system by 
the high frequency of near-equal sized binaries, and by 
the presence of binaries with extremely wide separations 
and long mutual-orbit periods. Widely-separated, long- 
period TNBs arc difficult to create and very sensitive to 
perturbation (Nesvorny et al. 2011, Parker & Kavelaars 
2010, Petit & Mousis 2004), and make valuable tracers of 
the dynamical and coUisional conditions over the history 
of the outer solar system. The orbital, compositional and 
statistical properties of these binaries constrain the total 
mass and dynamical history of the various populations, 
with important implications for theories of Solar System 
formation and evolution. 

There are several dozen known TNBs, but only a 
small subset has measured orbital parameters (Noll et 
al. 2008b, Naoz et al. 2010, Grundy et al. 2011). Most 
of these, in turn, are relatively tightly-bound binaries 
that have been characterized by observations from space 
(systems with mutual semi-major axis much less than 
5% of their HiU radius, eg., Grundy et al. 2009 & 2011). 
Two TNB systems with moderately widely-separated 
components have published mutual orbits (1998 WW31 
and Tcharonhiawako/Sawiskera, both with mutual semi- 
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major axes of order 5% of their Hill radius), but the 
widest TNBs (those with mutual semi-major axes sub- 
stantially exceeding 5% of their Hill sphere) have not 
been well-characterized to date, with a preliminary orbit 
estimate available only for the system 2001 QW322 (Pe- 
tit et al. 2008). Such wide-separation, near-equal mass 
binaries all have low heliocentric inclinations, indicating 
that they belong to the cold component of the Classi- 
cal Kuiper Belt. These wide binaries make up at least 
1.5% of the known Cold Classical belt objects (Lin et al. 
2010). 

The creation of large-separation, near-equal mass bina- 
ries is most likely during the formation phase of TNOs, as 
most proposed formation scenarios require a much higher 
space density of objects than is observed today. Ad- 
ditionally, most TNBs have identically-colored compo- 
nents, with differences between primary and secondary 
color being much smaller than the large color variation 
seen across the entire population of TNOs (Benecchi et 
al. 2009). This suggests that each component formed 
from similar material in a similar region of a locally ho- 
mogeneous protoplanetary disk with global variations in 
composition. Different mechanisms proposed for binary 
formation dominate under different dynamical conditions 
(eg., Schlichting & Sari 2008a). If the dynamical proper- 
ties of the systems today can be taken to be representa- 
tive of their primordial distribution, they can probe the 
dynamical conditions of the primordial Kuiper Belt dur- 
ing the formation phase. However, any intervening vio- 
lent dynamical events, such as collisions (Petit & Mousis 
2004, Nesvorny et al. 2011) or close encounters with gi- 
ant planets (Parker & Kavelaars 2010) can leave today's 
mutual orbit distribution substantially altered from its 
original state. It is critical to measure the orbital proper- 
ties of a large sample of TNBs, as well as perform dynam- 
ical studies of possible sources of orbital modification, in 
order to understand the full extent of information about 
the formation and history of the outer Solar System en- 
coded in these systems. 

We have collected astrometric measurements of a sam- 
ple of seven of the widest-known TNBs for an extended 
period, covering four to nine years of orbital motion for 
each system. These observations have allowed us to com- 
pute accurate mutual orbits for our sample of ultra-wide 
TNBs, and from these orbits we derive system mass and 
a host of other characteristics. In the first part of this 
paper, we outline the nomenclature we adopt to describe 
these systems and their host populations (§1.1), our sam- 
ple selection criteria (§2), details of our observational 
campaign and data reduction techniques (§3), and mu- 
tual orbit fitting algorithm (§4). In the latter part, we 
describe the mutual orbit fits (§5) and compare them to 
the properties of other binary populations, and derive 
geometric albedos for each system given reasonable as- 
sumptions of bulk density (§6). Finally, we conclude with 
a discussion of possible formation mechanisms and im- 
plications for the early history of the outer solar system, 
susceptibility of these systems to disruption by collisions 
and Neptune scattering, and present future surveys' abil- 
ities to discover and characterize a large sample of these 
ultra-wide TNBs. 



1.1. Nomenclature 

In this paper, we compare several sub-populations of 
Trans-Neptunian Objects and their various orbital prop- 
erties. In order to facilitate a clear understanding of the 
nomenclature we use to describe these populations and 
their properties, we provide an outline here. 

A binary's mutual orbit properties will be described ci- 
ther as a "mutual" property or denoted by the subscript 
"m." In contrast, the properties of the orbit of the bi- 
nary's barycenter around the Sun will be described as an 
"outer" property or denoted by the subscript "out." 

In order to compare the properties of our sample with 
those in the literature, some dynamical classification is 
required. We adhere roughly to the Gladman et al. 
(2008) nomenclature when discussing outer orbit proper- 
ties. In this paper, we frequently deal with binaries that 
belong to the following dynamical classes: 

• "Classical:" Non-resonant objects in the range 34 
AU< (7out < 47 AU, 37 AU < Cout ^ 70 AU. 

• "Cold Classical:" Subset of "Classical" objects 
with low orbital excitations and confined in semi- 
major axis. When dividing samples, we assign 
"Classical" binaries with Iq^^ < 10°, qQ■^^f^ > 38 
AU, and 42.4 < a^-^^f^ < 47 AU to this population. 
Referred to as CC population in the text. 

• "Hot Classical:" Subset of "Classical" objects with 
higher mean orbital excitations, and an extension 
to lower pericenter than the CC population. When 
dividing samples, we assign "Classical" binaries 
with iout > 10°' 9out < 38 AU, a^^t < 42.4 or 
"out > 47 AU to this population. Referred to as 
HC population in the text. 

This dynamical classification is somewhat different 
from that adopted by Grundy et al. (2011), and several 
binaries in that work which were classified as "extended 
scattered" fall into our HC classification. 

In addition, we compare the outer orbital distributions 
of binary sub-samples with the Canada-France Ecliptic 
Plane Survey (CFEPS) L7 synthetic model of the Kuiper 
BellQ. We compare the CC binary sub-sample with the 
composite of the "stirred" and "kernel" sub-components 
of the synthetic Kuiper Belt model, and refer to the com- 
posite of these sub-components as CC-L7. We compare 
the HC binary sub-sample with the "hot" sub-component 
of the synthetic Kuiper Belt model, and refer to this sub- 
component as HC-L7. 

In reality, any simple inclination cut is insufficient to 
determine which population a given object truly belongs 
to, as both the CC and HC populations overlap signif- 
icantly. According to the CFEPS L7 model, most rel- 
atively bright objects below 10° of inclination actually 
belong to the HC-L7 population. We stress that while 
we will refer to a given object as a "CC" binary or an 
"HC" binary, there is no way to absolutely verify the 
parent population for a given single object. However, 
we show later that the sub-samples of binaries which fall 
into our CC and HC classifications have dynamically dis- 
tinct mutual orbit distributions, and the outer orbit dis- 
tributions of CC binaries suggest that they are in fact 



^ Available at |http : //www. cfeps .net/L7Release .html | 
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members of the CC-L7 population, and likewise the HC 
binaries' outer orbit distribution is consistent with the 
HC-L7 distribution. 

2. SAMPLE SELECTION 

Since we seek to characterize the widest binaries (which 
have correspondingly long mutual periods), we opted to 
pursue a ground-based observation campaign. We chose 
our sample based on the following criteria: 

1. The system had no well-characterized orbit in lit- 
erature. 

2. The separation at discovery was > 0".5. 

3. The magnitude difference between the system's pri- 
mary and secondary was less than 1.7, indicating a 
near-equal mass system {Mp/AIg < 10). 

At the time of our sample selection, there were seven 
systems that met these criteria: 2000 CF105, 2001 
QW322, 2003 UN284 , 20 05 EO304, and three objects dis- 
covered over the course of CFEPS, with internal desig- 
nations b7Qa4, L5c02, and hEaV. The binary nature of 
2000 CF105 was presented in Noll et al. (2002), while 
the binary natures of 2003 UN284 and 2005 EO304 were 
presented in Kern (2006). Provisional orbital charac- 
terization for 2001 QW322 was presented in Petit et al. 
(2008). The CFEPS target L5c02 was identified as bi- 
nary by Lin et al. (2010), and the binary nature of b7Qa4 
and hEaV arc presented here for the first time. All of 
the outer orbits of this sample of objects fall into our CC 
classification, and have very low outer inclinations and 
eccentricities. 

Two other CFEPS targets, L4ql0 and L4kl2, were 
initially included in our sample, due to data from 
the Canada- France- Hawaii Telescope (CFHT) suggesting 
that they were elongated in a manner consistent with a 
near-equal mass binary with a separation of order 0".5. 
However, follow-up observations in very good seeing did 
not bear out their putative binary nature, and they were 
removed from our target list. The three CFEPS objects 
in our sample also have MPC designations: b7Qa4 is 
2006 BR284, hEaV is 2006 JZgi, and L5c02 is 2006 CHgg. 
Throughout this paper we will refer to these three sys- 
tems via their MPC titles. 

3. OBSERVATIONS AND DATA REDUCTION 

A targeted observational campaign from 2008 — 2011 
was executed from Gemini North using the Gemini 
Multi-Object Spectrograph in imaging mode, taken with 
the ^(50303 filter. Observations were queue scheduled, 
with stringent requirements on image quality (frequently 
at the expense of photometric conditions). By requiring 
modest visit times (~ 30 minutes), excellent seeing could 
be obtained without the use of adaptive optics, in some 
cases with full-width half-max of F ~ 0".35 or better. 
Additional observations during this period were made 
from VLT with the F0RS2 instrument, though image 
quality requirements were not held to the same stringent 
limits. Single-epoch observations were also made in April 
2010 from Magellan with the Megacam imager. 

Significant archival data also exists for all systems. We 
used the Solar System Object Search (S.D.J. Gwyn 2011, 



in prep) service provided by the Canadian Astronomy 
Data Centr^ to locate and download images from the 
CFHT and HST public archives that contained our tar- 
gets, and we also located images of our targets from the 
Mayall, Hale, and WIYN telescopes. 

Literature astrometric measurements are available for 
some systems as well. The relative astrometry for 2001 
QW322 published in Petit et al. (2008) is also included 
in our fit for that system. Astrometric measurements 
of 2003 UN284 and 2005 EO304 were presented in Kern 
(2006), and we include those measurements in our fits 
for these systems. 

Combining all these data sources, the smallest number 
of observed epochs for any binary in our sample is 12 vis- 
its for 2005 EO304, while the largest number is 35 visits 
for 2001 QW322- During most visits, more than one us- 
able image was acquired. The number of visits and total 
number of images from which astrometric measurements 
were made are listed in Table 1. 

Astrometric solutions were generated for each individ- 
ual image, matched to the J2000 coordinate system us- 
ing reference stars in the USNO b astrometric catalog. 
Whenever possible, the catalog stellar positions were cor- 
rected for proper motion since their last observed epoch, 
and uncertainties in the final reference positions reflected 
the original astrometric precision and the integrated un- 
certainty due to the stated uncertainty in proper motion 
over the intervening time. 

The brightest 100 non-saturated stars were identified in 
the CCD on which the binary was located (in the case of 
multi-chip imagers, other chips in the array were not used 
to constrain the astrometric solution), and their {x,y) 
positions (and uncertainties) extracted using SExtractor 
(Bertin & Arnouts, 1996). In the case of an image with 
a good initial World Coordinate System (WCS) estimate 
(eg., Gemini GMOS images), this WCS was used to es- 
timate each stars' RA and DEC position in the J2000 
system and the nearest neighbor in the USNO b catalog 
was identified as its matching counterpart. 

If an image did not have a good initial WCS, a robust 
pattern-recognition algorithm identified probable rough 
corrections to the WCS, applied these corrections, and 
then identified nearest-neighbor stars in the USNO b cat- 
alog. 

Once USNO b (RA, DEC) positions were matched to 
(x, y) positions in the image, the IRAF package ccmap 
was used to generate a WCS solution. Because this 
package does not handle positional uncertainties in ei- 
ther the image or reference positions, the positions of 
each matched star is cloned 1,000 times, adding Gaussian 
noise to each position consistent with the (RA, DEC) 
and (a;, y) uncertainties. Iterative fitting followed by au- 
tomatic clipping of outliers in these thousands of cloned 
sources allowed ccmap to generate a robust astrometric 
solution automatically which reflected the uncertainties 
in the absolute and measured positions of the reference 
stars. This allowed a more robust automatic solution to 
be derived with little input from the operator for each 
image processed. 

After the first pass of ccmap, more matches are 
searched for in the image with the USNO b catalog refer- 
ence stars, and upon fiagging any new matches the ccmap 

^ |http: //«ww2 . cadc-ccda.hla-lha.nrc-cnrc . gc . ca/ssos/| 
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Table 1 

Observations and System Properties 



Name 


Date Range 


^ visits 


Nobs 


IT^sys 


Am 


Hp'' 


Outer Orbit 












ir') 


ir') 


ir') 


«out (AU) 


^out 


W (°) 


2000 CFios 


2002-2011 


12 


50 


23.85" 


0.72(5) 


7.70 


43.84 


0.0362 


0.528 


2001 QW322 


2001-2010 


35 


88 


23.16'' 


0.03(5) 


7.51 


43.98 


0.0242 


4.808 


2003 UN284 


2003-2010 


14 


60 


22.7= 


0.88(6) 


7.5 


42.62 


0.0035 


3.069 


2005 EO304 


2005-2011 


12 


52 


22.45^* 


1.45(3) 


6.59 


45.62 


0.0679 


3.415 


2006 BR284 


2006-2011 


20 


66 


23.0° 


0.50(4) 


7.3 


43.80 


0.0393 


1.157 


2006 JZsi 


2006-2011 


15 


56 


22.7° 


0.98(2) 


6.9 


44.70 


0.0804 


3.550 


2006 CH69 


2004-2010 


15 


47 


23.0° 


0.44(5) 


7.0 


45.74 


0.0362 


1.791 



Notes. 

°: From CFHT MegaPrime, using Elixir photometric solutions. 
From Gemini North, Petit et al. (2008). 
From Gemini North, this work. 

From KPNO Mayall telescope, Benecchi et al. (2009). 
Assuming phase correction of 0.14 magnitudes per degree. 



Data Data-Model 




SEP: 1.59+/-0.06-, PA: 97.1" -10 -5 5 10 

FWHM 4.336 px, 0.632" 



Fig. 1. — Example of Gemini data and PSF fit. Top left: Original 
image from GMOS camera of 2001 QW322. Bottom left: Synthetic 
PSF model of binary components. Top right: Image residuals af- 
ter subtracting binary and other point-sources in the image. Bot- 
tom right: Relative contributions of both PSF components. Same 
stretch is applied to all images, and flux scaling is linear. 

routine is ealled again. This matching and WCS-fitting 
process is iterated ten times for every frame. 

The lowest-order astrometric solution merited by the 
distortions of the optics of each imager was used in each 
case. In the case of Gemini GMOS images, this was a 
simple rotation and fixed pixel scale. In the case of most 
other imagers, the distortion across the area of a single 
CCD was low enough such that the addition of indepen- 
dent X and y pixel scales, as well as a skew term, was 
sufficient. In the case of HST observations, the inter- 
nal HST astrometric solutions and distortion corrections 
were used. 

Once an astrometric solution had been found for an 
image, the relative astrometry of the binary pair in that 
image was then extracted using a custom Point-Spread 
Function (PSF) fitting routine. The PSF model we 
adopted was a sum of two elliptical Gaussian compo- 
nents, with arbitrary long-axis orientation for each com- 
ponent. The wider of the two Gaussians (the PSF wings) 
is arbitrarily limited to contain less than 2/3 the flux of 
the narrower "core" Gaussian to prevent runaway so- 
lutions with extremely wide wings. Adopting a non- 
circular PSF was required due to the fact that most im- 



ages were obtained with sidereal tracking rates, and over 
the period of integration the PSF of the binary's compo- 
nents became somewhat stretched along their direction 
of motion. 

A variant of the same algorithm used for fitting the 
mutual binary orbits (described in the following section) 
was used to minimize the residuals in a sky-subtracted 
40x40 pixel region centered on the binary. An initial 
interactive step is used to identify all the point sources 
in this region and flag the two associated with the bi- 
nary. Initial estimates for the amplitude and F of the 
PSF are made automatically, and these values are fed 
into a Markov- Chain Monte Carlo algorithm which finds 
the PSF model and array of point source positions and 
amplitudes that produces the smallest residuals. Be- 
cause the image model varies with the number of point 
sources in the 40x40 pixel region, the minimum num- 
ber of free parameters the algorithm must search over is 
9 (two point sources, and a two-component PSF model 
forced to be circular) while the maximum number of free 
parameters ever treated was 22 (five point sources and a 
two-component PSF model with arbitrary rotation and 
cllipticity for each component). Example of data from 
the Gemini observatory is illustrated in Figures 1 and 
2, along with the image residuals after PSF fitting and 
subtraction, and the model of the binary system. Each 
fit is visually inspected, and in general we found that our 
adopted PSF model produced extremely low residuals. 

In some cases where the two binary components were 
blended, we performed a check to verify that the ex- 
tra degrees of freedom added by allowing the PSF to be 
elongated was not skewing the measured astrometry. In 
these cases, we performed a second, independent fit using 
a circular PSF and compared the measured astrometry 
for the binary components. In general, we found excel- 
lent agreement between the two. In cases of very strong 
blending, fits with the elongated PSF would occasionally 
have trouble converging to a stable solution, and in these 
cases we adopted the values from the circular PSF fits. 
We did not attempt to determine any upper limits on 
separation based on completely unresolved images, and 
therefore at present such images do not contribute to our 
mutual orbit fits. 

Relative astrometry was recorded as separation (in arc- 
seconds) and position angle (in degrees East of North), 
and the observation date was taken at the central Julian 
Date of each observation. Uncertainty in the relative 
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astrometry was estimated as axy — F^SN^^ + SN^^, 

where SN indicates signal-to-noise ratio of primary or 
secondary. A noise floor is set by the uncertainty in the 
astrometric solution for each frame. 

These PSF fits also returned relative photometry for 
each system, and the mean A-mag measured in the Gem- 
ini GMOS rcosos filter during well-resolved visits is also 
included in Table 1. Relatively few observations were 
made in photometric conditions, as image quality was 
our primary concern. All MFC targets in our sample 
have published absolute system photometry in various 
bands, but only 2001 QW322 and 2005 E6304 have r'- 
band magnitudes in literature. 2000 CF105, 2006 BR284, 
2006 JZgi, and 2006 CHgg were all imaged on photo- 
metric nights from CFHT, and Elixir processed images 
were used to determine r'-band system magnitudes for 
these systems. The r'-band magnitude of 2003 UN284 
was determined from observations on a single night from 
Gemini North, though the absolute calibration of these 
particular images is poor and the resulting photometric 
uncertainty is relatively large. 

4. MUTUAL ORBIT DETERMINATION 

The basic operations performed by our mutual-orbit 
fitting routine are, given an initial guess of mutual orbit 
properties, to solve Kepler's equation in order to deter- 
mine the relative system geometry at the time of obser- 
vation (accounting for the light-time delay between the 
system and observer), then rotate the system in space to 
account for its orientation with respect to the Ecliptic. 
Finally, the code "observes" the system by applying a 
second rotation to account for the variation in viewing 
geometry induced by the relative motion between the 
Earth and the binary, and projects the result onto the 
sky plane given the separation between the observer and 
the system. 

To fit our observations, we chose to adopt the Metropo- 
lis algorithm minimization routine (Metropolis et al. 
1953), using a similar implementation to that described 
in Simard et al. (2002), who utilized the Metropolis algo- 
rithm to fit a 12 dimensional bulge -I- disk model to im- 
ages of galaxies. This algorithm is robust to complicated 
topology in parameter space, and can be easily adjusted 
to thoroughly explore parameter space at the expense 
of speed. The Metropolis algorithm is a Markov-Chain 



Table 2 

Mutual Astrometry 



Name 


JD 


Sep 


PA 


Err 




(d) 


(") 


{°) 


(") 


2001QW322 


2452117.92296 


3.6056 


290.0270 


0.2030 


2001QW322 


2452145.74085 


3.8915 


291.3190 


0.1722 


2001QW322 


2452145.76181 


3.8683 


287.1030 


0.1845 


2001QW322 


2452145.78289 


3.9600 


290.0990 


0.1168 


2001QW322 


2452146.75689 


3.7377 


292.2420 


0.2255 



Table 2 is published in its entirety in the electronic 
edition of the Astrophysical Journal. A portion is 
shown here for guidance regarding its form and content. 



Monte Carlo technique which, after an initial burn-in pe- 
riod, occasionally makes "bad" decisions, allowing it to 
diffuse out of local minima. After a number of iterations, 
the choice of new parameter values can be informed by 
previous values, improving the speed of convergence in 
complicated parameter-space topology. A binary mutual 
orbit has seven free parameters, and in our implemen- 
tation wc chose these to be the following: mutual semi- 
major axis {am), eccentricity (em), period (T^), Mean 
Anomaly (M, valid at a defined JD), inclination (i^), 
longitude of the ascending node {^e), and argument of 
pericenter (we). The last three angular parameters are 
defined with respect to the J2000 Ecliptic. For nearly 
circular orbits, M and u) become degenerate and alter- 
nate choice of basis is preferred; however, we were not 
presented with a circumstance where altering the basis 
used in our code was merited, as the binaries we ob- 
served all have significant eccentricity. For all orbit fits, 
15 Metropolis algorithm threads are run simultaneously, 
and each compare their final best-fit value and sampled 
orbit space to identify the global best-fit and statistically 
acceptable range of parameters after an additional test 
of the error distribution. 

Because our estimates of the astrometric uncertainty of 
each observation may not reflect their true distribution, 
(i.e., we don't know the properties of the distribution of 
errors on the individual measurements a priori), we per- 
form a test to bootstrap the 68% and 95% confidence in- 
tervals for the distribution. This test simply assumes 
that the best-fit orbit is the "truth," and estimates the 
observed error distribution around the best-fit. The test 
randomly draws n measurements of the observed error 
from the best-fit orbit from the pool of n real observa- 
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tions of a given target (sampled with replacement, so 
observations may be repeated in the resampled list). Af- 
ter this resampling, we compute a new "observed" chi^ 
based off of the resampled list and best- fit orbit. We 
store this new chi^ statistic and repeat the process 10,000 
times, building up a distribution of best-fit chi^ values 
which "might have been." We find the 68% and 95% 
upper confidence intervals on this distribution, and set 
those to be our thresholds for statistically acceptable 
mutual orbit fits. Orbits that fall below these thresh- 
olds are saved, and their distribution is used to generate 
the uncertainties for each orbital parameter. This anal- 
ysis has led to extremely well-behaved fitting behavior, 
with consistently nested uncertainty contours after every 
addition of new data to a given fit. 

The mutual orbit fitting code was tested by reproduc- 
ing the orbital parameters for the Pluto-Charon system 
based on synthetic data generated by the JPL Horizons 
system, and reproducing the orbital parameters of 2001 
XR254 and 2004 PB204 as published in Grundy et al. 
(2009) to an accuracy well within their stated uncertain- 
ties. We did not explore the robustness of our method 
for determining confidence intervals with these tests, as 
we did not have sufficient numbers of real observations 
for the two published mutual orbits we tested. 

5. PRESENT BEST-FIT ORBITS AND IMPLICATIONS 

The astrometric measurements for each system and 
comparison to fit orbits are illustrated in Figures |3] and 
m These figures project each system onto the sky plane 
in physical separation units, removing the variation in 
the observed scale of the systems due to the change in 
the observer-system separation over the course of a year. 
An example of the mutual astrometry for all the TNBs 
used in these Figures is shown in Table 2. The full table 
is available in the electronic edition. 

The best-fit and uncertainty (given as the extrema of 
each parameter from the distribution of orbits allowed 
at the 68% level of confidence) of all fitted mutual or- 
bit parameters are listed in Table 3. Additionally, Table 
4 contains derived parameters; specifically, the system 
mass Msys, mutual semi- major axis to Hill-radius frac- 
tion a/Rn, mutual inclination im, mutual argument of 
pericenter ujm, and mutual pericenter separation in mul- 
tiples of primary radii Qm/Rp- Figure [5] illustrates the 
a/Rn, e™, and im fits and the 68% and 95% uncertain- 
ties in these parameters for each system. 

5.1. Derived Parameters 

In the following section, we describe the derived mutual 
orbital parameters listed in Table 4 and illustrated in 
Figure[5] System mass is simply calculated from Kepler's 
laws, 

M.ys=^n^^. (1) 

while the classical Hill radius for a binary in orbit around 
the Sun is defined as 

1 

Rh = aout(l - eout) [jm^) ' 
where Oo^t ^'^'^ ^oui '^^ ^^'^ heliocentric semi-major 



axis and eccentricity of the binary system's barycenter. 
Primary radius is found from the system mass by as- 
suming that both components have the same albedo and 
density, and is given by 

Rp = f ^„ -3An. ,^'-^SyS ) , (3) 

y47rp(l + 10^^) ) 

where p is the adopted bulk density and Am is the magni- 
tude difference between the primary and secondary com- 
ponents of the binary. 
Mutual inclination is the angle between the pole vector 

of the binary mutual orbit Pm and that of the outer orbit 

Pouti ^1^*^ '^^^ ^6 found by i,„ = cos^^(Pout ' Pm)- The 
pole vectors for cither orbit can be found by (eg., Naoz 
et al. 2010): 

(sin(ri) sin(i) \ 
- cos(17) sin(i) . 
cos(i) ) 

The mutual argument of pericenter Um (critical for es- 
timating the extent of Kozai oscillations) is the angle 
between the ascending node (with respect to the outer 
orbit) and pericenter in the plane of the mutual orbit. It 
can be found by: 

ujra cos"^ ( ,1"' , , " ) sign(em [z] ) , 

where 

^ / cos(u;_E)cos(ri£;)-sin(w£;)cos(i_E)sin(f2E) \ 
-rr-,= cos(iU£;)sin(f2_E)+siii(w_E)cos(z£;)cos(f2£;) 
1^"! V sin(wi5) sin(iB) j 

and n = Pout Pm- 

5.2. Kozai Cycles 

Systems within a broad range of im centered on 90° 
may be subject to large oscillations in im and due 
to the Kozai effect (Kozai 1962). Over the oscillation 
cycles two values are conserved: one depending on the 
initial mutual eccentricity and inclination, and the other 
depending on both these and the mutual argument of 
pericenter. Following Perets & Naoz (2008), we adopt 
the following form for these two conserved values: 

A = (5e^, sin^ ujm + 2(1 - e^)) sin^ i„ (4) 

and 

B = ^y^-el^ cosim- (5) 

With some algebraic manipulation and the constraint 
that eccentricity and inclination minima and maxima oc- 
cur when u!m = 0° or 90°, these two constants determine 
the maximum and minimum eccentricity and inclination 
reached by a given binary system during its Kozai cy- 
cle. We calculate these values, and the amplitudes of the 
eccentricity excursions experienced by each binary sys- 
tem is listed in Table 5. Also listed are the predicted 
minimum pericenter separations (occurring during the 
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Fig. 3.— Astromctry and fitted mutual orbits for MPC binaries 2000 CF105, 2001 QW322, 2003 UN284, and 2005 EO304. Latitude 
separation (A0) and longitude separation (AA) are given in projected physical units (km) to remove variation due to changing separation 
between binary system and the observer and illustrate physical scale of each system. Black line indicates best-fit mutual orbit, while dark 
and light gray regions illustrate orbits consistent at the 68% and 95% confidence level, respectively. Mutual astrometry is available online 
as a machine-readable table. 



Table 3 

Fit Mutual Orbit Elements 



M 

n 

262. 

158+11 

294+I4, 

193.1?3 

219.ti 

104.19 

286.±^ 



Name 



(lO"! km) 



(years) 



n 



n 



2000 CF105 
2001 QW322 

2003 UN284 

2005 EO304 

2006 BR284 
2006 JZgi 
2006 CH69 



0.06 
+0.38 
0.14 
t; cc + O.SS 
'^•'^'^-0.53 

6.98t«i? 

r, c-q + 0.03 

^•■^■^-0.03 
o i-}Q-|-0 . 53 
■^•^■'-0.28 
1-0.33 
0.28 



3.33 
10.15 



2.76^ 



10.92 
17.01 
8.73 



0.10 
1.55 

0.69 
+0.65 
0.54 
+0.45 
-0.45 
+0.04 
0.03 
-,-1+0.15 
'-'--0.12 
oQ+0.05 
°^-0.07 



ir 



0.29 
0.46 
0.40 
0.22 
0.28 
0.84; 
0.90 



TTTUT- 

0.02 
+0.02 
0.01 
+0.04 
0.07 
+0.02 
0.02 
+0.01 
0.01 
+0.03 
0.02 
+0.02 
0.02 



*~0.7 

150.7«j 

24.3^1-5 
1+0.8 
-0.5 
55.6ljj 
13.3+^' 



167.4"' 



12.4"* 



134.1 



1.9 
+4.9 
-6.1 



223. 
243 

92 
259 

41. 

82 
105. 



3 
+3 
-4 
+6 
-3 
+2 
3 
2 
2 
+5 
-7 
+6 
8 



n 



Epoch 



296.+3 

172.1^0 
206.^5 
14.ti 
171. 
149.t^ 



2454880.96 
2452117.92 
2452963.77 
2453440.94 
2455153.08 
2455007.86 
2455190.06 



highest eccentricity phase of the Kozai cycle) in multi- 
ples of primary radii, as close encounters during Kozai 
cycles may lead to modification of the mutual orbit due 
to tidal friction (Fabrycky & Tremaine 2007, Perets & 
Naoz 2009, Brown et al. 2010). Details of systems with 
large- amplitude Kozai oscillations will be discussed on 
an object-by-object basis in the following sections. 

The Kozai effect can be easily suppressed by other 
effects, including permanent asymmetries in the mass 
distribution of the component bodies (Ragozzine 2009). 
However, this suppression only occurs for relatively small 
semi-major axes, and therefore the Kozai oscillations of 
these wide systems are unlikely to be suppressed. 



5.3. Individual Objects 
5.3.1. 2000 CF105 

With several epochs of HST data, and nearly nine years 
of observational baseline (covering most of a single 11 
year mutual orbit period), 2000 CF105 has one of the 
best-measured orbits in our sample with mutual semi- 
major axis and period uncertainties at the 1% level. It 
is also the lowest-mass system in our sample, and the 
second-most weakly bound (behind only 2001 QW322)- 
Because of its low mass (currently the lowest mass of 
any known TNB), this system has the smallest estimated 
primary radius of any system in our sample (assuming 
all objects share a common density), estimated to be 
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Fig. 4.— Same as Figurc[3l but for the three CFEPS binaries 2006 BR284 {b7Qa4), 2006 JZgi (hEaV), and 2006 CHeg (L5c02). Mutual 
astrometry is available online as a machine-readable table. 



Table 4 

Derived Values 



Name 


-^'^s y S 

(1017 kg) 


a-m/RH 


(°) 


n 


qm/Rp* 


2000 CF105 

2001 QW322 
2003 UN284 

2005 EO304 

2006 BR284 
2006 JZgi 
2006 CH69 


J-°^-0.14 

21.50_2 23 
1 q 1 0+2.26 

_L0.1Z_2 gy 

21.03t°-?I 

'^-0.20 
8.30l:^;35 


1679+"'""1^ 

n QOOO~l~^ ' ^ ■^'^^ 

U.ZZZZ_Q QQgJ 

n 1449+°°°™ 

U.1449_g QQgQ 

n 1 c:c:q+0.004S 

U.1000_Q QQ47 

U-U8'9_Q 0005 

0900+'' °°21 
u.uyuu_Q QQjg 

0.0809_o 0009 


167.9+^!^ 
22.7+1-4 

i5.7i«;« 

54.6lli 

11 1+2-5 
^^■^-2.0 

133.3ltl 


295. +3 

248.i?o 
165.121 
203.+g 
13.+1 
158.1* 
147.11 


74l+*J 

aKc:+64 
o00_40 

534lf2 

7141^ 

408+? 

84l}°s 
56ll« 



Primary radius R,p assumes p = 1 gram cm 



31.8+Q-g km given a bulk density of 1 gram cm"''. 

2000 CF105 is not subject to strong Kozai oscillations, 
and its pericenter separation is al-ways greater than ~ 710 
primary radii, so there is little chance of mutual tides 
having modified its mutual orbit. 

The pole solution for 2000 CF105 is non-degenerate at 
greater than 95% confidence, and the system is retro- 
grade. Its mutual pole vector is only ^ 12.1° degrees 
from anti-aligned -with its outer orbit's pole vector, mak- 
ing it one of the most pole-parallel systems kno-wn. 

5.3.2. 2001 QWi22 

At the outset, our current results for the mutual orbit 
of 2001 QW322 appear inconsistent -with the orbit pub- 
lished in Petit et al. (2008). By fitting only the astromet- 



ric data used in that paper "with our ne-w code, we find 
a much larger range of allo-wable mutual orbit solutions 
than was presented in the previous study. This is espe- 
cially notable for allowing much higher values of mutual 
eccentricity. This larger range of allowed orbital param- 
eters completely overlaps with our current orbit fit, as 
illustrated in Figure El We suspect that the difference is 
due to two factors: our more thorough fitting algorithm 
and our different statistical analysis of the allowed 
range. The allowed range in Petit et al. (2008) was 
much smaller than that determined by the methods used 
in this work. 

The components of this system remain photometrically 
indistinguishable, with Am consistent with 0. When 
making astrometric measurements, we have arbitrarily 
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Fig. 5.— Best-fit mutual orbit properties. Blue: 2000 CFips, Red: 2001 QW322, Yellow: 2003 UN284, Purple: 2005 EO304, Green: 2006 
BR284i Gray: 2006 JZgi, Brown: 2006 CHgg. Heavy and light contours represent 68% and 95% confidence intervals, respectively, while 
black stars mark the best-fit parameters. Teal symbols mark sample of synthetic systems created by gravitational collapse in Nesvorny 
et al. (2010), selecting simulations with Q = 0.5 — l.OQcirc (initial clump rotation), and /* = 10 — 30 (cross-section modifier). Gray 
shaded region represents properties of binary systems created by three-body exchange reactions in Funato ct al. (2004). Green shaded 
region represents properties of binary systems created by Chaos- Assisted Capture in Astakhov et al. (2005). Other points represent TNBs 
with known orbits (Grundy et al. 2011). Black points are members of the Classical belt, gray points are members of other dynamical 
populations. Filled points indicate Am < 1.7, and open points indicate Am > 1.7. Upward triangles indicate CC objects, while downward 
triangles indicate HC objects. Square points indicate resonant objects, and diamond points indicate Centaurs or scattered disk objects. 
Non-preferred degenerate pole solutions are not illustrated for literature binaries. 



assigned the Northernmost object in the discovery epoch 
as the system primary. Observations have been frequent 
enough that there is no possibility for confusion between 
system components, based on continuity of the orbital 
motion. 

The observations in our dataset cover approximately 
nine years, sampling slightly over 50% of the best-fit mu- 
tual orbital period of 17.0 years. Because the observa- 
tions have been frequent and of high quality, this limited 
sample of the orbital motion of the system is very con- 
straining. Mutual semi-major axis, period, and eccen- 
tricity are all known to better than 10% accuracy. An- 
gular and derived parameters are similarly well known. 

2001 QW322 is subject to strong Kozai osciUations, 
with the nominal best-fit orbit implying eccentricity vari- 
ations between 0.342 < e™ < 0.477. Its current orbit is 



therefore near its highest eccentricity phase. 

In physical units, this system remains the most widely 
separated binary minor planet known, with a mutual 
semi-major axis of 1.015lQ g^4 x 10^ km. It is also likely 
the most weakly bound binary minor planet known, with 
its measured a/i?jy of 0.2222i°;oJgi exceeding the current 
estimate for the outer satellite of the Main-Belt aster- 
oid (3749) Balam of a/Rn ~ 0.2 (Marchis et al. 2008). 
Several other known main-belt asteroid binaries have es- 
timates of a/Rn shown in Richardson & Walsh (2006) 
which are similar to or slightly higher than that mea- 
sured for 2001 QW322, but these estimates are based on 
single-epoch observations and may not refiect the true 
orbits and masses of these systems. 

The pole solution for 2001 QW322 is non-degenerate 
at greater than 95% confidence, and the system is retro- 
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Fig. 6. — Comparison of the allowed orbits for 2001 QW322, 
given the data presented in Petit et al. 2008 (gray contours) and 
the complete set of data used in this work (black contours) using 
the fitting algorithm described in the text. Stars mark best-fit 
orbits. 68% and 95% confidence intervals are marked by the solid 
and dashed contours, respectively. 

grade with a mutual inclination of ^ 152.7°. 

5.3.3. 2003 UN2S4 

The orbit of 2003 UN284 is the least-well constrained 
in our sample. The fit relies heavily on astrometry pub- 
lished in Kern (2006) for pinning the 2003—2004 astrom- 
etry, and only three data points constrain the 2005 — 2008 
astrometry. Recent data has proved relatively discrimi- 
natory, and the mutual semi-major axis and period are 
both known to better than 10% accuracy, but the de- 
rived system mass is only constrained to ~ 20% accuracy. 
The Am we adopt for 2003 UN284 is determined from 
two well-resolved visits from Gemini North. Kern (2006) 
finds a highly variable Am for this system, suggesting 
one or both components may have significant lightcurve, 
and our adopted Am may not reflect this variability. 

This system is likely subject to minor Kozai oscilla- 
tions, but the amplitude of these oscillations are compa- 
rable to the uncertainty in the currently measured mu- 
tual eccentricity. The pole solution for 2003 UN284 is 
non-degenerate at greater than 95% confidence, and the 
system is progradc with a mutual inclination of ^ 22.7°. 

5.3.4. 2005 EOzm 

The orbit of 2005 EO304 also relies heavily on astrom- 
etry published in Kern (2006), and it has the fewest ob- 
served epochs (12) of any of our systems. Nevertheless, 
the recent measurements from VLT and Gemini have 
provided reasonably tight constraints on the orbit prop- 
erties. 

This system has the largest Am in our sample at 1.45 
in r' . It also has the largest primary radius (assuming 
all objects share a common density), which we estimate 
to be 76.21q g km given a bulk density of 1 gram cm~^. 

2005 EO304 is not subject to significant Kozai oscilla- 
tions, and its minimum pericenter separation is at least 



692 primary radii. The pole solution for 2005 EO304 is 
non-degenerate at greater than 95% confidence, and the 
system is prograde with a mutual inclination of ^ 15.7°, 
making it one of the lowest mutual inclination systems 
known. 

5.3.5. 2006 BR2Si 

Nearly all observations of 2006 BR284 used in our as- 
trometric fit come from Gemini^ and the orbit has been 
sampled for just over a single orbital period. Its orbit is 
very well constrained, with mutual semi-major axis and 
period known to 1% accuracy and mutual eccentricity 
to better than 4% accuracy. The pole solution for 2006 
BR284 is non-degenerate at greater than 95% confidence, 
and the system is prograde with a mutual inclination of 
~ 54.6°, making it the most inclined system in our sam- 
ple. 

2006 BR284 is subject to strong Kozai oscillations, with 
the nominal best-fit orbit implying eccentricity variations 
between 0.72 < < 0.26. Its current orbit is therefore 
near the lowest point of its eccentricity cycle. Its pericen- 
ter passages are still widely separated (never lower than 
143 primary radii) and mutual tides are not a concern 
for this system. 

5.3.6. 2006 CHeg 

The observations of 2006 CHeg cover well over a single 
orbital period, though as Figure S] illustrates its on-sky 
behavior is quite complex and well time-sampled observa- 
tions were necessary to accurately constrain this system's 
mutual orbital properties. 

This system has the highest well-measured mutual ec- 
centricity of any binary minor planet, at e„i = 0.901q'q2- 
The orbit of the outermost satellite of the trinary as- 
teroid (3749) Balam has been estimated to rival this at 
Gm ^ 0.9 (Marchis et al. 2008), but its value is poorly 
constrained. A fascinating consequence of this extreme 
mutual eccentricity is that over the course of its mutual 
orbit, the secondary of 2006 CHgg will subtend an angle 
ranging from to ~ 0.1° (at mutual apocenter) to ~ 1.7° 
(at mutual pericenter) as viewed from the surface of the 
primary — from one-fifth to over three times the angular 
size of the Moon on the sky. 

2006 CHgg may be subject to strong Kozai oscillations, 
with the nominal best-fit orbit implying eccentricity vari- 
ations between 0.70 < e„i < 0.94. Its current orbit is 
therefore near its highest eccentricity phase. The pole 
solution for 2006 CHgg is non-degenerate at greater than 
95% confidence, and the system is retrograde with a mu- 
tual inclination of ~ 134°. This system is also the most 
tightly bound in our sample, with a/Rn — 0.0809. Be- 
cause of the orbits' high eccentricity, the error distribu- 
tions in am and Tm conspire to produce a large uncer- 
tainty in the derived system mass, and Adgys remains 
uncertain at the 40% level. 

The high mutual eccentricity phases of 2006 CHgg's 
Kozai cycles lead to very close passages between the pri- 
mary and secondary. At its current mutual eccentricity, 
pericenter passages occur at 56^]^g primary radii (again 
given a bulk density of 1 gram cm~'^), which may be wide 
enough that mutual tides do not cause orbital modifica- 
tion. However, during the high eccentricity phase of its 
Kozai cycle, the pericenter separation of 2006 CHeg drops 
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Fig. 7. — Features of 2006 CHsg's Kozai oscillations. Top panel: 
representative sample of 107 orbits consistent with 2006 CHgg as- 
trometry at Icr level, showing their current eccentricity (eo) and 
the maximum eccentricity (emax) they reach over the course of 
a Kozai cycle, with the best-fit orbit marked by the large point. 
Bottom panel: Same sample of orbits, but now illustrating their 
best-fit system mass (Mgyg) and the minimum pericenter separa- 
tion in multiples of primary radii (qmin/ Rp, assuming p = 1 gram 
cm~'^). 

much lower to Sl^y. Figure [7] illustrates the distribution 
ol mutual pericenter separation versus system mass. If 
we argue that 2006 CHgg must have survived roughly in 
its current orbital configuration for the age of the solar 
system, then we would prefer higher minimum pericenter 
separations to keep the system from experiencing tidal 
evolution. From Figure [7] we see that orbit fits which 
have higher pericenter separations are lower-mass solu- 
tions. However, the tidal evolution of highly eccentric, 
highly inclined binary systems is poorly understood, and 
future work is merited to determine if limiting the tidal 
evolution of 2006 CHgg would provide useful priors for 
further constraining its mutual orbit. 

5.3.7. 2006 JZsi 

The observations of 2006 JZgi span slightly more than 
a single period, and Figure |4] illustrates that like 2006 
CHgg the on-sky behavior of this system is also com- 
plex. This system is also very highly eccentric, at 
= 0.84 

-01)2' making it the second- highest eccentricity 
TNB known behind 2006 CHgg. Its mutual semi-major 
axis remains somewhat poorly constrained at 16% un- 
certainty, while the mutual period is known to better 



than 5% uncertainty. The derived system mass remains 
highly uncertain for the same reason as 2006 CHgg, and 
uncertainties in Mgys remain at the 40% level. 

The pole solution of 2006 JZsi is non-degenerate at 
the 95% level, and the system is prograde with the low- 
est mutual inclination of any known TNB at just ^--^ 11°. 
Due to this low inclination, 2006 JZgi is not subject to 
notable Kozai oscillations. Despite its high eccentricity 
and relatively large primary (61^g° km), mutual pericen- 
ter passages are always at least 64 primary radii, making 
this system much less susceptible to possible mutual tidal 
effects than 2006 CHgg. 



5.4.1. 



5.4. Ensemble Results 
Comparison to Other Populations 



Comparing the TNBs studied here to previously char- 
acterized TNBs, we see that they are much more widely 
separated in terms of their Hill sphere occupation. Fig- 
urcOillustrates the mutual semi-major axis to Hill radius 
fraction for the ultra-wide TNBs characterized in this 
work and the systems listed in Grundy et al. (2011); 
the widest known TNB with a well-characterized or- 
bit not in our sample is Teharonhiawako/Sawiskera at 
q/Rh ^ 0.06. 

Grundy et al. (2011) showed that the previously 
known widely-separated, loosely-bound TNBs arc only 
found on dynamically cold heliocentric orbits. Here we 
seek to confirm this relationship, and identify which dy- 
namical population hosts the wide binaries; is it the cold 
Classical Kuiper Belt, or can a low-inclination exten- 
sion of the hot Classical Kuiper Belt plausibly host the 
widely-separated binary systems? Figure[S]illustratcs the 



outer orbital excitation (given by y sin(iQy^)2 + e^^^) 

versus a/Rn- As in Grundy et al. (2011), we con- 
firm that only dynamically cold populations host wide 
binaries. Further, systems with pericenters which sug- 
gest current or past encounters with Neptune have rel- 
atively low a/Ru, supporting the destructive nature of 
such encounters presented by Parker & Kavelaars (2010). 
To quantify the difference between the binaries found in 
dynamically cold and hot populations, we compare the 
a/Rn distributions of binaries falling into our CC clas- 
sification with all other binaries. We find that these two 
samples are inconsistent with being drawn from the same 
a/Rn distribution, with the KS test rejecting this hy- 
pothesis at greater than 99.9% confidence. Therefore, the 
binaries which fall into our crude dynamical classification 
of the cold Classical Kuiper Belt have distinctly different 
characteristics than binaries hosted in other populations. 

To further clarify the host population of the wide bina- 
ries, we compare the outer orbital excitation distribution 
of the binaries with the distribution of the CC-L7 and 
HC-L7 orbital excitations. The distribution of orbital 
excitations for these binary samples and the L7 model 
populations are illustrated in Figure [9] Since most bina- 
ries considered here have been discovered in near-ecliptic 
surveys, there is significant bias against the detection of 
high inclination binaries. The bias is less significant when 
considering low inclination populations, so comparisons 
between the biased CC binary sample with the de-biased 
CC-L7 model population are reasonably fair. Comparing 
the CC binary sample to just the low inclination objects 
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Fig. 8. — Left panel: Heliocentric orbital excitation vs. a/Rn (similar to Fig. 5 from Grundy et al. 2011). Triangle, square, and 
diamond points have same meaning as in Figure [S] representing orbits published in Grundy et al. (2011). Stars mark best-fit orbits for the 
ultra-wide TNBs characterized in this work, with crrorbars representing 68% confidence interval. Vertical lines mark upper 95th-percentile 
of the orbital excitation of the cold Classical Kuipcr Belt (solid) and hot Classical Kuiper Belt (dashed), as found by the CFEPS L7 survey. 
Right panel: heliocentric pcriccntcr qout vs. a/Rjj. 



in the de-biased HC-L7 model population is also mean- 
ingful, as observational biases are not significant for the 
low-inclination extension of the HC-L7 population. We 
use the KS test to determine if we can reject the fol- 
lowing two hypotheses: that the CC binaries are drawn 
from the CC-L7 orbital excitation distribution, or that 
the CC binaries are drawn from the excitation distribu- 
tion of low-inclination (^out — 5°) HC-L7 orbits. Addi- 
tionally, we determine whether just the seven ultra-wide 
binaries characterized in this work (with no outer orbit 
constraints placed on them) could be drawn from the 
CC-L7 or Iq^^ < 5° HC-L7 excitation distributions, and 
whether all binaries with a/Rn > 0.02 and Am < 1.7 
could be drawn from the same excitation distributions. 

We rule out (P < 0.001) that the CC binaries are 
drawn from any subset of the HC-L7 excitation distri- 
bution, while it is plausible (KS test cannot reject at 
high confidence) that they are drawn from the CC-L7 
excitation distribution. We also find that, without any 
prior cuts on heliocentric orbits, neither the seven ultra- 
wide binaries characterized in this work nor all binaries 
with a/ Rh > 0.02 and Am < 1.7 can be drawn from the 
low-inclination HC-L7 excitation distribution (P < 0.001 
in both cases), while both can be plausibly drawn from 
the CC-L7 excitation distribution. This further confirms 
that wide binaries are intimately linked to the dynami- 
cally cold population of the Classical Kuiper Belt — it 
cannot be that they are predominantly hosted by a low- 
inclination extension of the hot Classical Kuiper Belt. 

5.4.2. Mutual mclmation distribution 

The mutual inclination of a TNB system is one of the 
most challenging parameters to measure, as there is a 
mirror degeneracy in the pole solution which can only 



Table 5 

Kozai Oscillations 



Name 


emax 






2001 QW322 
2003 UN284 
2006 BR284 
2006 CHgg 


'-0.006 

'-'■^3-0.07 
n, 79 + 0.02 
'^-0.02 
94+0 01 
^■="^-0.02 


n qq+0.04 
26-^+° °°^ 

U.ZOO_g Qj^Q 

70+0 06 
^•"^-0.08 


830+29 
4861^^ 
155111 
31+? 


2000 CF105 

2005 EO304 

2006 JZgi 


^■^*-0.03 
*J-°ti-0.02 


n 90+0.02 

U-^°-0.02 
n. 99+0.02 
'-'■^■'-0.03 

n 040.03 

U-»4_o 02 


7391^9 
698+g 
82llg 



* : Primary radius Rp assumes p = 1 gram cm . 

be broken after sufficient time has elapsed for the the 
observer's viewing geometry of the binary system has 
changed enough to discern the system's true orientation. 
However, the distribution of mutual inclinations and the 
ratio of prograde to retrograde orbits holds significant 
implications for formation scenarios (Sehlichting & Sari 
2008b, Noll et al. 2008b) and for the ongoing evolu- 
tion of the binary orbit (Fabrycky & Tremaine 2007, 
Perets & Naoz 2009). As illustrated in Figure El aU 
seven systems characterized in this work now have non- 
degenerate pole solutions: four prograde and three ret- 
rograde. The prograde-to-retrograde ratio and its 95% 
Poisson counting uncertainty for the ultra-wide TNBs 
is therefore ^ 1.33^^' 02. If we fold in the systems with 
non-degenerate pole solutions presented in this work and 
those in Grundy et al. (2011) which fall into our CC 
sub-sample and meet our near-equal mass criteria, we 
find that the ensemble prograde-to-retrograde ratio for 
dynamically cold, near-equal mass TNBs is l.GOtlil. 

The distribution of the mutual orbit inclinations of the 
ultra-wide binaries is illustrated in Figure [TOl along with 
the inclination distribution of tightly-bound TNBs in lit- 
erature, taken from the orbits compiled by Grundy et al. 
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Fig. 9. — Cumulative histograms of heliocentric orbital excitation for different binary sub-samples, compared to model distributions. 
Left panel: All CC binaries (heavy solid histogram), just the ultra-wide binaries (gray histogram), and all binaries with a/Rjj > 0.02 and 
Am < 1.7 (light black histogram), compared to the CC-L7 model distribution (dashed line). Right panel: All HC binaries (solid histogram) 
compared to the HC-L7 distribution (dashed histogram). We stress that this plot is comparing a heavily biased sample (binaries) with a 
de-biased model, and is only used for consistency checks. 




Fig. 10. — Cumulative mutual inclination distribution of ultra- 
wide binaries characterized in this work (heavy black histogram), 
tight binaries from literature (light black histogram), and tight 
CC binaries from literature (dotted black histogram). Retrograde 
inclinations are folded over onto the range of 0° — 90°. Dashed 
black line shows ultra-wide binaries' inclinations when smoothed 
over entire Kozai cycle of each object. Gray dashed line: uniform 
distribution — P(i) oc sin{i). Grey dash-dotted line: uniform 
distribution modified by removing any orbits with Kozai cycles 
which drive q^^ij^ < 30Rp as described in the text. Gray solid line: 

inclination distribution of form P(i) oc sin(i)e~''-^('/^^ ) . 



(2011). Due to the relative lack of high- inclination and 
excess of low-inclination wide binaries, we find it highly 
unlikely that the wide binaries' inclinations are drawn 
from the same inclination distribution as the tighter 
TNBs, with the Kolmogorov-Smirnov (KS) test reject- 
ing this hypothesis at > 95% confidence. Additionally, 
we find that the tight binaries in literature are consis- 
tent with having their inclinations drawn from a uni- 
form distribution (P(i) oc sin{i), with the KS test re- 
jecting this hypothesis at only ^ 6% confidence), while 
the ultra-wide TNBs' inclinations are inconsistent with 
being drawn from the same uniform distribution (KS test 
rejecting this hypothesis at ~ 99% confidence). A sim- 
ple way to understand the strength of this rejection is 
to note that in a uniform inclination distribution, 50% 
of the inclinations arc > 60°, while there are no systems 
in our sample with inclinations so high. Since the prob- 
ability of randomly drawing an inclination < 60° is 0.5 
each time, after sampling seven systems the probability 
of every system having inclination < 60° is the same as 
landing seven coin flips head-up in a row, 0.5^ ~ 0.008. 

When restricting the tightly-bound literature binaries 
to only those in the CC sample, the significance of the 
rejection of the wide binaries being drawn from the same 
distribution drops to only > 89% — however, the dy- 
namically cold binaries from literature remain consistent 
with a uniform distribution (KS test rejecting this hy- 
pothesis at only 40% confidence). 

It should be noted that observational bias works 
against the detection of low inclination systems (with re- 
spect to the Ecliptic), as they present edge-on geometry 
and their average projected separations are lower than 
high inclination systems. Thus, our relative lack of high 
inclination objects is not due to observational bias. 
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When considering all the inclinations the wide binaries 
reach during their Kozai cycle, we find that a sin(i) times 
a Gaussian distribution (frequently used to describe the 
inclination distribution of heliocentric Kuipcr Belt orbits, 
eg.. Brown 2001) centered at i = 0° with width a ~ 22° is 
compcUingly similar to the Kozai-modificd distribution. 
We find that widths between 10° < cr < 50° are consis- 
tent with the observed (non-Kozai) distribution within 
the 95% confidence interval. 

The difference in inclination distributions between 
both populations could either be due to cosmogonic vari- 
ations inherent in different formation mechanisms, or due 
to evolutionary processes. Fabrycky & Tremaine (2007) 
showed that tidal friction coupled with Kozai oscilla- 
tions can cause wide, high-inclination systems to shrink 
and become circularized, creating a paucity of widely- 
separated, high-inclination binaries, and an excess of 
tight binaries near the critical inclinations for Kozai cy- 
cles (40° and 140°). This is referred to as the Kozai 
Cycles with Tidal Friction (KCTF) mechanism, and its 
plausibility as a significant evolutionary mechanism for 
binary minor planets was confirmed by Ragozzine (2009) 
and Perets & Naoz (2009). While small numbers of ob- 
jects and other complicating effects may prohibit the 
detection of an increase of tight binaries near the criti- 
cal inclinations, the observed paucity of high-inclination, 
ultra-wide TNBs is suggestive of the KCTF effect in ac- 
tion. 

We have performed a cursory test to compare the ob- 
served inclination distribution to the outcomes of KCTF. 
We determine the maximum eccentricity Cmax each bi- 
nary system can reach for a grid of initial inclinations 
io and arguments of pericenter wq, assuming they be- 
gin with initial e™ and equal to their present values. 
We then determine the fraction of phase space as a func- 
tion of iQ (assuming iq was initially uniformly distributed 
on a sphere) which do not lead to an emax which cause 
each system's pericenter to drop below a critical num- 
ber of primary radii — in other words, we require that 
Qmin > n X Rp. Ragozzine (2009) showed in full numer- 
ical simulations with reasonable assumptions that tidal 
dissipation became significant at q ^ 20Rp for the bi- 
nary system Orcus/Vanth, which is more massive than 
the binaries studied here. However, the binary system 
2006 CHgg has Kozai cycles which take it to 31~t^Rp, 
and presumably its existence suggests that such pericen- 
ter separations are stable for a significant fraction of the 
age of the solar system. As such, we chose the limit 
qmin > 30 X Rp for our cursory KCTF test, and compare 
the resulting "KCTF-modified" inclination distributions 
to the observed ultra-wide TNB inclination distribution 
(illustrated in Figure [TU)) . We find that the resulting dis- 
tribution (and any with smaller qmin cutoff) is ruled out 
at 97% confidence. 

We conclude that while the KCTF mechanism may 
have modified the orbits of some very high-inclination 
ultra-wide TNBs, it alone is not enough to explain the 
current mutual inclination distribution of these systems. 
To verify this, more comprehensive studies of the effects 
of KCTF on systems like those presented here are needed. 
If KCTF is not sufficient to explain the mutual inclina- 
tion distribution of the ultra-wide TNBs, then we must 
consider cosmogonic effects. The primordial poles of the 
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FlG. 11. — Albedos and radii for CC binary systems. Left/bottom 
axes show values assuming p = 1 gram cm 
show the values assuming p = 0.6 gram cm 
mary radii, triangles connected with dashed line mark secondary 
radii. Black points represent binaries characterized in this work, 
gray points represent values derived from literature for the low- 
inclination classical TNBs 1998 WW31, Borasisi/Pabu, Logos/Zoe 
and Teharonhiawako/Sawiskera. Dotted lines represent contours 
of constant H magnitude. Gray region marks H,. > 8, likely un- 
populated due to flux limits of current binary searches. 

wide binaries may have preferred orientations orthogo- 
nal to the Ecliptic plane, suggesting formation in a very 
cold disk (Noll et al. 2008b). It should be noted that 
if the inclination distribution of the ultra-wide TNBs is 
non-uniform, and the binaries are subjected to small per- 
turbations over their lifetimes (e.g., collisions), then the 
mutual inclination distribution will always tend to be- 
come more random over time, approaching a uniform 
distribution. As such, the primordial inclination distri- 
bution would have to have favored low inclinations even 
more strongly than the current distribution does. We 
explore the effects of collisions on the mutual inclination 
distribution in an upcoming paper. 

6. ALBEDOS AND DENSITIES 

Given the dynamically-derived masses and visible pho- 
tometry for each system, we can explore the albedos 
and densities for the component bodies in these bina- 
ries. Without radiometric measurements to ascertain 
the albedo independently, the density and albedo remain 
degenerate. However, by assuming physically plausible 
values for the component densities, we can estimate the 
implied albedos for each object; the results of this ex- 
ercise are illustrated in Figure [11] and listed in Table 6. 
Generally, albedos for our sample of ultra-wide TNBs are 
found to be consistent with those measured radiometri- 
cally for larger solitary cold Classical Kuiper Belt objects 
(eg., Brucker et al. 2009), and range from 9% — 30%, as- 
suming (0=1 gram cm~^ (6.4% — 21% assuming p = 0.6 
gram cm~^). Figure [Til also includes estimates of the 
albedos and radii of four other binaries from literature, 
and these systems were selected as members of the CC 
sample which had estimates of their r-band magnitude. 
The albedos of these four literature systems range from 
5.4% — 28% (with p=l gram cm~^). 

The apparent trend of increasing albedo with decreas- 
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Table 6 

Albedos and Primary Radii 
(with p = 1 gram cm~^) 





p 




Note 




(r') 


(km) 




9000 PPinr 


'-'■■^'-'-0.03 


31 8+"-'^ 
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1 


2003 UN284 
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'°-^-0.9 




2006 BR284 




44.91°;^ 




2006 JZgi 


17+0-07 


6111° 




2006 CHeg 


n 90+0.09 

^■^'-'-0.06 


50+1 




1998 WW31 




741^3 


1,3,7 


Teharonhiawako 


n 147+0.003 

'J-J^^' -0.003 


76.0l°j 


1,4,7 


Borasisi 


„ ,00+0.009 

U.iyZ_g QQg 


81.6l°i 


2,5,7 


Logos 


n 90+0.02 


41.lt«i 


2,6,7 



1; Adopting m© = -26.93. ^: Adopting mo = -27.12. ^: 
Photometry from Veillet et al. (2002). ^: Photometry from 
Benecchi et al. (2009). ^: Photometry from Delsanti et al. 
(2001). ^: Photometry from Jewitt & Luu (2001). ^: Mass from 
Grundy et al. (2011). 

ing radius visible in Figure [TT] is due to a selection effect: 
in a flux-limited survey, any physically small object de- 
tected must have a high albedo. The apparent trend here 
is consistent "with a flux limit some"what less than Hr — 8. 
We note that, in general, the observations which discov- 
ered the primary of a given binary system "were not the 
observations "which discovered the binary nature of the 
system. Thus, this flux limit seems to apply to the pri- 
mary absolute magnitude, and the fact that secondary 
absolute magnitudes scatter across the Hr — 8 line re- 
flects deeper follo"w-up observations identifying the sec- 
ondaries. 

The lack of a strong group of small, high-albedo bi- 
naries along "with the lack of large, high-albedo binaries 
suggest that relatively lo"w albedos may be more common 
than high albedos in the cold Classical Kuiper Belt. All 
seven binaries in Figure [11] "with primary radius greater 
than 55 km have nominal albedos less than 20%, "while 
all four binaries "with primary radius smaller than 55 km 
have nominal albedos greater than 20%. Given the steep- 
ness of the size distribution of low-inclination TNOs in 
this size range (g ^ 4.8, Fraser & Kavelaars 2009), we 
would expect there to be roughly nine times as many ob- 
jects in the radius radius range from 30 — 55 km than at 
all radii larger than 55 km. Since the sample contains 
only roughly 0.57 times as many binaries in the smaller 
size range, we posit that the rest of the expected small 
binaries are missing due to having low albedos, making 
them invisible to the flux limits of the surveys which dis- 
covered the binary systems in our sample. 

We adopt an ansatz albedo distribution of a Gaussian 
centered a.t p = 0.05, and clipped such that p > 0.05 
(comparable to the lowest measured albedo in our sam- 
pie): 

{f. r-^ P-Q.05 -,2 

Additionally, we adopt a size distribution with slope 
q = 4.8, and estimate the flux limit at discovery to be 
Hr = 8. When we draw a radius from this size distri- 
bution, we assign it an albedo from our Gaussian albedo 



distribution and determine if it is bright enough to have 
been observed by our synthetic survey (brighter than 
Hr = 8). We then compare the properties of these syn- 
thetic "observed" systems with the systems which were 
actually observed; we vary the width Up of the albedo 
distribution until the KS test can rule out that either 
the distribution of real radii or the distribution of real 
albedos are drawn from the synthetic "observed" distri- 
butions. We find that at 95% confidence, the observed 
range of real albedos implies that the albedo distribution 
must have a width ap > 0.058 or there would be too few 
high albedo detections. Additionally, given the observed 
real primary radius distribution, the albedo distribution 
must have a width Up 0.1 or there would be too many 
small objects discovered. Figure [T2] illustrates the distri- 
bution of albedos and primary radii of the real binaries 
in comparison to the synthetic "observed" distributions 
of these parameters. 

The lack of a strong group of small, high-albedo bina- 
ries may also be explained in the case of a more uniform 
albedo distribution by positing that the binary fraction 
decreases drastically for decreasing radii, similar to the 
prediction of Nesvorny et al. (2011). However, with the 
addition of a varying binary fraction with size as a new 
degree of freedom, the current sample size is not suffi- 
cient to quantitatively constrain such behavior at this 
time. We note that given the large range of albedos ob- 
served in this population, any sharp features in the trend 
of binary fraction with radius will be blurred if consid- 
ering only the absolute magnitude of systems (as done 
in Nesvorny et al. 2011), and any such features will be 
much more evident when radii derived from mass mea- 
surements or radiometric measurements are used in place 
of absolute magnitudes. 

7. DISCUSSION 

7.1. Formation Mechanisms and Implications 

Since the discovery of the flrst TNBs, a number of pos- 
sible formation pathways have been posited. In the fol- 
lowing discussion, we consider those mechanisms most 
likely to form widely-separated, near-equal mass systems, 
and compare the predicted outcomes of each of these 
pathways to our observed sample. 

7.1.1. 1/2 s and L3 mechanisms 

Originally described in Goldreich et al. (2002), these 
mechanisms were further investigated by Noll et al. 
(2008b) and Schhchting & Sari (2008 a & b). The LjS 
mechanism posits that binaries are captured when two 
passing solitary objects can disperse some excess kinetic 
energy into a sea of smaller bodies and become bound, 
while the L3 mechanism instead sends the excess kinetic 
energy away through scattering a third large body. 

Schhchting & Sari (2008b) show that models like L2S, 
which rely on a smooth dissipation process to capture 
binaries, will dominate the binary formation rate only 
in conditions where the relative velocity between plan- 
etesimals v is much less than the Hill velocity, vh = 
2ttRh /Tq^^, where Tq^^^ is the heliocentric orbital pe- 
riod. They also show that under these conditions, the bi- 
nary mutual inclinations will be dominantly retrograde, 
predicting a prograde-to-retrogradc ratio < 0.03. The 
measured wide binary inclinations exclude such an ex- 
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Fig. 12. — Left panel: comparison of binary albedo distribution (heavy histogram) with model distributions, assuming an albedo distri- 
bution with form of Eqn. (6] a power-law size distribution with slope q = 4.8, and binary discovery survey flux limit of Hr = 8. Intrinsic 
albedo distributions shown by dashed lines, while synthetic "observed" distributions shown by solid lines. Right panel: Same binaries and 
albedo distributions, but showing resulting radius distributions. 



treme ratio of prograde to retrograde systems. There- 
fore, we can rule out this mechanism for forming wide 
binaries, unless an intervening dynamical process can 
be invoked to rc-orient a large number of binary sys- 
tems or preferentially destroy retrograde binary systems. 
We estimate that starting from a primordial prograde- 
to-retrograde ratio of 0.03, at least 22% of wide binary 
systems would have to be re-oriented in order to not be 
ruled out at greater than 95% confidence by the current 
observed prograde-to-retrogradc ratio. 

Under more energetic conditions, where the relative 
velocity between planetesimals exceeds vh, Schlichting 
& Sari (2008b) show that three-body interactions (L3 
models) will dominate the binary formation rate. In this 
regime, they find that roughly equal numbers of pro- 
grade and retrograde systems are formed, consistent with 
the observed distribution of inclinations. However, they 
also show that only systems with separations of order 
s < RnivH/vY tend to survive the formation phase, 
and binary formation rates drop dramatically as v in- 
creases. As our observed binary systems have separa- 
tions of order 0.08 — 0.23i?H! the velocity dispersion in 
the primordial disk could not have exceeded 2 — 4 times 
vh if this formation mechanism applied. Additionally, 
Noll et al. (2008b) suggest that formation in a dynam- 
ically cold disk {v < vh) should produce aligned orbit 
poles. Since the wide binaries seem to prefer low mutual 
inclinations, this argues for formation in a dynamically 
cold disk, but not so cold as to allow L2S to dominate 
and produce a large fraction of retrograde systems. 

Together, the widely-separated components, lack of 
clear preference for retrograde orbits, but apparent pref- 
erence for low mutual inclinations all point toward the 
velocity dispersion being approximately equal to the Hill 
velocity. This represents a fine-tuning problem (Noll et 



al. 2008b), for there is no clear a priori reason to expect 
that V ~ Vh- Additionally, it is not clear whether the 
balance between the L3 and L2S mechanisms at w ^ vh 
would simultaneously produce widely separated binaries, 
aligned poles, and roughly equal numbers of prograde 
and retrograde orbits. 

7.1.2. Exchange Reactions and Chaos-Assisted Capture 

Funato et al. (2004) suggest that multiple exchange re- 
actions (where one object in a binary system is swapped 
for a passerby) can produce very widely separated bi- 
nary systems. However, systems as widely separated as 
those in our sample formed through exchange reactions 
all have very high eccentricities (e™ > 0.9, see Figure[5]). 
The systems 2006 CHeg and 2006 JZgi have present ec- 
centricity values consistent with these predictions, but all 
other systems are presently inconsistent with such high 
eccentricities. Several systems are subject to large oscil- 
lations of their inclination and eccentricity due to Kozai 
cycles, but these oscillations do not carry them to ec- 
centricities as high as predicted by exchange reactions. 
Thus, it seems unlikely that this mechanism dominated 
binary formation. 

Astakhov et al. (2005) simulated the effect of chaotic 
transient binaries on stable binary formation. They 
found that two objects temporarily caught in their mu- 
tual chaotic layer could become stabilized by dynami- 
cal friction due to a sea of small objects — effectively 
adding an enhancement to the L2S mechanism due to 
transient, chaotic orbits. This mechanism is referred to 
as Chaos- Assisted Capture. They find mutual eccentric- 
ities spanning the range of those observed in our sam- 
ple, but the separations they find for binaries formed by 
Chaos- Assisted Capture do not extend to as high as those 
found for the ultra- wide TNBs (see Figure [5l). Addition- 
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Fig. 13. — Cumulative histogram of a/Rn for CC TNBs (solid 
histogram) and all other TNBs with known orbits (dashed his- 
togram). Also shown is results of gravitational collapse binary 
formation simulations by Nesvorny et al. (2010), selecting simula- 
tions with n = 0.5 — l.Oficirc (initial clump rotation), /* = 10 — 30 
(cross-section modifier), and final mass ratio < 10. 

ally, Schlichting & Sari (2008a) argue that the enhance- 
ment due to these transient captures is not significant, 
and that formation should proceed as they found for the 
L2S and I/3 mechanisms. We conclude that it is unlikely 
that this mechanism dominated the ultra-wide TNB for- 
mation rate. 

7.1.3. Gravitational collapse 

Recently, another mechanism has been proposed to 
form Kuipcr Belt binaries. Operating with the frame- 
work of planetesimal formation through rapid gravita- 
tional collapse in a turbulent disk, Nesvorny et al. (2010) 
suggest that binaries may form as a cloud of cm-scale 
particles collapses and fragments. This model produces 
binaries very efficiently, and their properties can vary 
widely. Mass ratios of order unity are produced, and 
semi-major axes from 10"^ — 10^ km are produced for sys- 
tems with primary radii ranging from tens to hundreds of 
kilometers. Broad ranges of inclination and eccentricity 
can be produced for all semi-major axes. Additionally, 
this mechanism has the attractive feature of producing 
a natural explanation for the correlated colors of binary 
components (Benecchi et al. 2009), in contrast to the 
broad range of colors exhibited between different binary 
systems. 

We compare the results of a subset of these simula- 
tions to our observed mutual orbits. We select only those 
simulations which produce binaries with final mass ra- 
tio < 10, and with initial particle-swarm rotation Q, = 

0.5 — l.Or^circ (where Vlcirc — \J GM , and M and 
R are the total initial mass and radius of the swarm) 
and collisional cross-section enhancements (to account 
for the lower resolution of the simulation compared to 



reality) /* = 10 - 30. See Nesvorny et al. (2010) for a 
more thorough description of these parameters and their 
importance. Figure [5] illustrates the mutual eccentric- 
ity and a/Rn for all the synthetic orbits formed by this 
mechanism which meets these criteria, and in general we 
find that they appear to mimic the distribution of real 
orbits surprisingly well. Figure [T3] illustrates the distri- 
bution of a/Rn for the synthetic orbits compared to the 
CC binaries and all the other binaries. We find that the 
a/Rn distribution of the CC binaries is statistically in- 
distinguishable from the synthetic orbit distribution gen- 
erated by gravitational collapse, while we rule out that 
the other binary populations were drawn from the same 
distribution as the synthetic orbits at a high level of con- 
fidence. 

This agreement is encouraging, but it remains prelimi- 
nary. The initial conditions for binary formation by grav- 
itational collapse is still highly uncertain, and it remains 
unclear if this mechanism can produce a large number 
of retrograde binaries or if the binaries it produces will 
exhibit a preference for low mutual inclinations. Addi- 
tionally, while the a/Rn distributions appear to agree, 
there are likely a number of biases in the current observed 
distribution and any conclusions drawn from it must be 
taken with caution. The binary fraction is measured to 
be ^ 30% in the cold Classical Kuiper Belt (Noll et al. 
2008a) at the limit of HST resolution, while the fraction 
of ultra- wide systems is approximately 1.5% (Lin et al. 
2010); therefore, the ultra- wide binaries make up approx- 
imately 5% of the current CC binary population. How- 
ever, comparing the systems with measured orbits, we see 
that the ultra-wide systems are over-represented, with 
our sample alone making up ~ 39% of all the CC binaries 
with measured mutual orbits. Taking this into account, 
it appears that the Nesvorny et al. (2010) model actu- 
ally over-produces wide binaries compared to the present 
population. 

There is additional concern since the current a/Rn dis- 
tribution may not represent the primordial distribution; 
Petit & Mousis (2004) showed that these wide binary 
systems were likely much more common in the past, and 
may have been reduced in number due to collisional dis- 
ruption. We discuss this problem in more detail in the 
following section. 

7.2. Susceptibility to Disruption 

Because of their low binding energy, ultra-wide TNBs 
are very susceptible to disruption by a variety of pertur- 
bations. Parker & Kavelaars (2010) simulated the trans- 
plantation of the Cold Classical Kuiper Belt from closer 
to the Sun via Neptune scattering as suggested by Levi- 
son et al. (2008), and found that wide binaries are very 
efficiently destroyed by close encounters with Neptune 
before they can be implanted in the current Kuiper Belt. 
The binaries we characterize in this work are all wide 
enough to be easily stripped by this process, and com- 
paring the a/Rn and eccentricity of our systems to the 
results of Parker & Kavelaars (2010), we find destruction 
probabilities ranging from at least 75% for 2006 BR284 
to over 98% for 2001 QW322- Furthermore, such interac- 
tions would likely have randomized the orbit poles of the 
surviving wide binaries, not leaving behind the aligned 
poles we see today. 

Taken together, we estimate that if our sample rep- 
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resents a population which has been subjected to dis- 
ruption by Neptune scattering, the initial population of 
ultra-wide TNBs would have to have been roughly 13 
times larger than the current population. If we take the 
lower limit of the current wide binary fraction of the Cold 
Classical Kuiper Belt estimated by Lin et al. (2010) at 
1.5%, this indicates that the primordial wide binary frac- 
tion have to exceed 20% to leave enough wide binaries 
surviving post-Neptune scattering. Correcting the esti- 
mate of Kern & Elliot (2006) to apply to only the Cold 
Classical Kuiper Belt (using the same approach as Lin 
et al. 2010) results in an estimate closer to 5% for the 
current wide binary fraction, which would imply a pri- 
mordial fraction of wide binaries in excess of 65%. These 
estimates assume no other processes disrupted binaries in 
the intervening time between implantation in the Kuiper 
Belt and the present day. A primordial wide binary frac- 
tion of 20% is comparable to the entire current binary 
fraction in the Cold Classical Kuiper Belt. 

As mentioned previously, the Nesvorny et al. (2010) 
simulations appear to over-produce wide binaries with 
respect to today's population; synthetic binaries as wide 
or wider than the ultra-wide TNBs characterized in this 
work created by the Nesvorny et al. (2010) model rep- 
resent roughly 20% of the systems produced by their 
simulations. Therefore, if the primordial binary fraction 
was ~ 100%, the ultra-wide binary population produced 
in these simulations would be roughly sufficient to leave 
a remnant population similar to the observed popula- 
tion after implantation in the cold Classical Kuiper Belt 
through Neptune scattering if no additional processes 
disrupted these systems afterward. However, since a 
much larger fraction of the tighter binaries would survive 
the implantation process, we would expect a significantly 
larger remnant population of tight binaries if the primor- 
dial binary fraction was so high. Later processes (such 
as collisions) could modify the binary fractions over the 
age of the solar system, but generally wide binaries are 
more susceptible to stripping processes and the relative 
fraction of wide to tight binaries will tend to decrease 
over time. Therefore, barring a formation mechanism 
which much more strongly favors wide separations, we 
find that the ultra-wide binary orbits confirm the find- 
ings of Parker & Kavelaars (2010) that it is unlikely for 
the cold Classical Kuiper Belt to have experienced a pe- 
riod of scattering encounters with Neptune and likely 
formed in situ. 

In addition to being susceptible to disruption by tidal 
effects, collisions are effective at unbinding these systems. 
Petit & Mousis (2004) estimated the mean lifetimes for 
wide binaries given an estimate of the number of im- 
pactors in the current Kuiper Belt capable of unbinding 
these systems in an impact. They found that binaries 
like 2001 QW322 had mean lifetimes much less than the 
age of the solar system, and required a primordial pop- 
ulation of wide binaries at least seven times larger than 
the extant population to explain the continued existence 
of several such systems. 

Following Petit & Mousis (2004) and A.H. Parker & 
Kavelaars (2011, in submitted), we estimate the mean 
lifetime of a binary on a circular mutual orbit due to 
coUisional unbinding by an impact on the secondary to 
be: 
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where A^o is the number of impactors larger than i?o, 
drawn from a power-law size distribution with slope q, 
Vi is the mean impact velocity, and Pi is the intrinsic 
coUisional probability. For the following discussion, we 
adopt the same values as Petit & Mousis (2004) for Pi = 
1.3 X 10^^^ km^^ yr^^, and Vi ~ 1 km s^^, which were 
based upon Farinella et al. (2000). 

This calculation can be somewhat refined by noting 
that collisions onto the primary can also contribute to 
unbinding the system, and when folding in this second 
decay channel the system's mean lifetime becomes half 
the harmonic mean of the individual coUisional lifetimes 
of the primary and the secondary, 



V,, 
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For equal- mass systems like 2001 QW322, this effec- 
tively halves the estimated mean lifetime. 

For an impactor population, we take the estimated to- 
tal population of TNOs in the CFEPS L7 model which 
cross through the cold Classical Kuiper Belt region and 
extrapolate to small sizes. To do this, we first estimate 
the number of objects at Hg = 10 (roughly the magni- 
tude of the putative break in the luminosity function of 
the Kuiper Belt) by extrapolating from the limits of the 
CFEPS survey {Hg ^ 8.5) with a luminosity function 
power-law slope of 0.76 (e.g.. Eraser & Kavelaars 2009). 
We estimate that there are of order Nq ^ 800,000 ob- 
jects brighter than Hg = 10 in which pass through the 
cold Classical Kuiper Belt. Translating this to a radius 
with a geometric albedo of 0.1, we adopt Rq ~ 26 km 
as our break and normalizing radius for the small-object 
size distribution. 

A number of groups (Bernstein et al. (2004), Eraser 
& Kavelaars (2009), Fuentes et al. (2009)) have found 
that the luminosity function of the Kuiper Belt breaks 
to a shallower slope of approximately a ~ 0.2, trans- 
lating to a power-law slope of q ~ 2. However, Eraser 
(2009) showed that the size distribution likely steepens 
again after a "divot," and the true small object size dis- 
tribution slope may be somewhat steeper than q = 2. To 
capture this uncertainty in the behavior of the size distri- 
bution below the break, we estimate the coUisional life- 
times of the ultra-wide TNBs using both a shallow slope 
of g = 2 and a coUisional-equilibrium slope of q = 3.5. 
This steeper slope is comparable to the estimate of the 
size distribution slope at small sizes set by the putative 
detection of a single sub-km TNO by stellar occultation 
presented in Schlichting et al. (2009). 

CoUisional lifetimes from Eqn. [5] for each binary sys- 
tem are presented in Table 7. This table also contains 
an estimate of the required initial number of binaries n 
to leave behind one system with a given set of properties 
after t = 4 x 10^ years of coUisional bombardment, given 
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Table 7 

CoUisional Lifetimes 



Name 


q = 2 




q = 3.5 






T (yr) 


n* 


r (yr) 


n* 


2000 CFio5 


5.0x10^" 


1 


1.0x10^ 


52 


2001 QW322 


2.6x10^" 


1 


2.7x10^ 


4 


2003 UN284 


3.3xlOio 


1 


2.5x10'' 


5 


2005 EO304 


3.1x101'' 


1 


2.9x10^ 


4 


2006 BR284 


4.2x10^'' 


1 


2.2x10^ 


6 


2006 JZgi 


3.7x101'' 


1 


3.0x10^ 


4 


2006 CH69 


3.9x101" 


1 


2.7x10^ 


5 



* : Estimate of primordial number of wide binaries with same 
system properties in order to leave one system surviving today, 
from n = e'*/'^) with t = 4 X 10^ years. 

n = e'*/"^-*. For the shallow small-object size distribu- 
tion slope g = 2, all the binary systems characterized in 
this work have collisional lifetimes much in excess of the 
age of the solar system, and if this slope represents the 
size distribution then the binary population we see today 
can be taken as representative of the primordial popula- 
tion. However, in the case of the steeper slope q = 3.5, 
we see that all the systems have lifetimes less than the 
age of the solar system, and each binary system repre- 
sents a member of a decayed initial population ranging 
from 4 — 52 times larger than the current population (the 
average required primordial wide binary population for 
this size distribution is ~ 11 times larger than the extant 
population). 

We conclude that if the power-law size distribution 
slope at small sizes remains shallow, then the extant wide 
binary population can be taken as representative of the 
primordial wide binary population. If this is the case, 
then the preliminary agreement between the predictions 
of binary formation by gravitational collapse and the ob- 
served orbital distribution of wide TNBs remains valid 
when considering the primordial wide binary orbital dis- 
tribution. 

The estimates of collisional lifetimes presented here arc 
very simple analytical estimates, and ignore some impor- 
tant effects such as the eccentricity of the binary system, 
mass loss during impacts, and solar tides. We have re- 
fined these calculations and performed extensive numer- 
ical simulations of binary disruption by collisions, and 
will present the results of this work in an upcoming pub- 
lication. 

The system 2000 CF105 has by far the shortest colli- 
sional lifetime, due to its wide separation and very small 
component sizes; in fact, it is the second-widest known 
TNB (with respect to a/Rn) and its components have 
the lowest mass of any TNO currently measured. It can 
be unbound by the smallest impactors of all the binary 
systems considered here, and is therefore most sensitive 
to the impactor population at small sizes. The frequency 
of systems like 2000 CF105 in the current Kuiper Belt is 
therefore of considerable interest, as they act as power- 
ful probes of the collisional environment. Future surveys 
will likely identify many such systems, and we discuss 
the prospects for these surveys in the next section. 

7.3. Processes in the primordial disk 

In the previous sections, we have found that the incli- 
nation distribution of the wide binaries indicates a very 
cold dynamical environment at the time of formation, 
but formation of the binaries post-accretion through the 



L2S pathway — which should dominate in such condi- 
tions — can be ruled out due to a lack of preference 
for retrograde orientations. We also found that the cold 
Classical Kuiper Belt was never subjected to a period 
of close encounters with Neptune, and therefore likely 
formed in situ. Finally, the properties of mutual orbits 
produced by current simulations of gravitational collapse 
are consistent with the observed properties of binaries in 
the cold Classical Kuiper Belt. 

Together, all these results suggest that formation of 
the cold Classical Kuiper Belt occurred in situ, with ob- 
jects forming via rapid collapse of small particles directly 
into large (10 — 100 km) planetesimals through primary 
accretion driven by turbulent concentration of solids in 
the protoplanetary disk (eg., Johansen et al. 2007, Cuzzi 
ct al. 2010). This process offers a mechanism for form- 
ing binaries such as are observed in this population, and 
can feasibly form the cold Classical Kuiper Belt in situ 
(Cuzzi et al. 2010) without requiring that a large fraction 
of the primordial mass in planetesimals be subsequently 
removed through either significant collisional grinding or 
dynamical processes (both of which can be destructive 
to the binary population), even when the density of the 
primordial disk is required to drop sharply beyond 30 
AU to halt the early migration of Neptune at its present 
distance. Additionally, since this process produces no 
massive sea of small planetesimals (most planetesimals 
being "born big"), the L2S mechanism will naturally be 
ineffective, thus preventing it from forming a dominantly 
retrograde binary population in the dynamically cold en- 
vironment indicated by the binaries' low mutual inclina- 
tions. 

The mechanism of binary formation by gravitational 
collapse should be extensively explored. If it can be 
shown to produce roughly equal numbers of prograde 
and retrograde binaries while maintaining a preference 
for low mutual inclinations, then this should be taken 
as strong support for the formation of the cold Classical 
Kuiper Belt in situ by rapid primary accretion driven by 
turbulent concentration and gravitational collapse. 

7.4. Characterizing wide binaries with next- generation 

surveys 

While the current sample of known ultra-wide TNBs 
remains fairly small, future wide-area surveys like the 
Large Synoptic Survey Telescope (LSST) and Pan- 
STARRS will be quite capable of identifying such sys- 
tems. The projected LSST g-band single-visit depth is 
^ 25, which translates to Hg ~ 8.5 for a geocentric dis- 
tance of 45 AU. At this depth, there are roughly 17,000 
objects in the CC-L7 component of the CFEPS synthetic 
model of the Kuiper Belt. At a minimum, the results of 
Lin et al. (2010) suggest that at least 1.5% of these ob- 
jects are ultra-wide TNBs with characteristics like those 
in our sample, and we therefore expect over 250 such 
systems will be detectable by the LSST pipeline. 

Some fraction of the observations of each of these TNBs 
will likely be resolved, depending on the binary separa- 
tion and the seeing at the time of each observation. In 
order to estimate this fraction of resolved observations, 
we have used clones of each binary in our sample to de- 
termine the distribution of on-sky separations that can 
be expected over the ten-year baseline of the LSST sur- 
vey. We move each clone to a random starting point on 
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Fraction of observations which resolve binary 

Fig. 14. — Results of observing clones of our sample of ultra- 
wide TNBs through simulated LSST seeing, assuming observations 
in the r-band taken over the course of the 10-year LSST survey 
lifetime. An observation is said to be "resolved" when the full- 
width at half-max of the seeing during that observations is smaller 
than the projected binary separation. 

both its heliocentric orbit and mutual orbits, then prop- 
agated its motion forward for ten years while sampling 
its on-sky separation at 100 "observation epochs" over 
than ten-year period. At every "observation epoch" we 
determined if the system was resolved or not based on 
its on-sky separation and an image quality drawn from 
the measured r-band distribution of seeing at the LSST 
site (LSST Science Book Version 2.0, Section 2.2), with 
median seeing of '--^ 0".6. Figure [Ml illustrates the cumu- 
lative fraction fs of binary systems versus the minimum 
fraction of observations which resolve those systems. 

The minimum fraction of resolved observations for any 
system we simulated was ~18%, and the maximum was 
100%. If we take the planned baseline number of 230 
visits per field in the r— band to represent the number 
of observations likely to be made of each binary, then 
we would expect ^40 — 230 resolved observations of each 
binary system. This indicates that virtually all 250 TNB 
systems with properties similar to those in our sample 
likely to be observed by the LSST will have more resolved 
observation epochs than the most-observed binary in our 
current data set (2001 QW322, with 35 observed epochs). 
As such, not only will LSST detect most of these binary 
systems, but it will measure their mutual orbit properties 
at least as well as we have managed to do with the sample 
presented in this paper. 

With such a large catalog of well-characterized binary 
orbits, all drawn from the same well-calibrated survey, 
unprecedented avenues of investigation will be opened. 
Identifying trends of binary properties with heliocentric 
orbit may provide further evidence for different origins 
and histories of the sub-components of the Kuiper Belt, 
and will provide new constraints on the dynamical his- 
tory of the outer Solar System (Parker & Kavelaars 2010, 



Murray-Clay & Schlichting 2011), trends of binary frac- 
tion with component sizes will constrain the extent of 
collisional grinding in the early Kuiper Belt (Nesvorny 
et al. 2011), and color — albedo trends may be identified 
and used to constrain the surface composition of these 
objects. 

8. SUMMARY 

We have presented the first-ever well-characterized 
mutual orbits for a sample of seven ultra-wide Trans- 
Neptunian Binaries. These orbits range over eccentric- 
ities from 0.2 — 0.9, and have semi-major axis to Hill 
Radius {a/Rn) fractions ranging from 0.08 — 0.22. We 
find that their properties are distinct from other binary 
populations, with the following highlights: 

1. The outer orbits of all widely-separated binaries 
{a/ Rh > 0.02) with near-equal mass components 
(Am < 1.7, or Mp/Ms < 10) are distributed con- 
sistently with being drawn from the "stirred" and 
"kernel" components of the CFEPS L7 synthetic 
model of the Kuiper Belt, and inconsistent with 
being drawn from the low-inclination subset of the 
"hot" component. This confirms that only the dy- 
namically cold components of the Classical Kuiper 
Belt are host to wide binaries. 

2. Two ultra- wide TNBs have very high mutual ec- 
centricities {cm — 0.84 and 0.9) consistent with the 
predictions of binary formation by exchange reac- 
tions (Funato et al. 2004), but the rest arc incon- 
sistent with the extreme eccentricities predicted by 
this mechanism. 

3. Ultra- wide TNBs have roughly equal numbers of 
prograde and retrograde orbits, and the observed 
ratio is inconsistent with the predicted preference 
for retrograde orbits if the L2S binary formation 
mechanism dominated unless over ^ 22% of wide 
binary systems have switched orientations over 
their lifetimes. 

4. Ultra-wide TNBs have a statistically significantly 
different mutual inclination distribution compared 
to tighter binaries in literature, and their inclina- 
tions also cannot have been drawn from a uniform 
sphere ( P{i) cx sin{i) ) due to a lack of detec- 
tions of high-mutual inclination wide binaries. This 
paucity of wide binaries at high inclinations and 
preference for low inclinations cannot be explained 
by the KCTF mechanism alone, and suggests a pri- 
mordial preference for low-inclination mutual or- 
bits. This suggests formation in a very dynamically 
cold disk. 

5. The wide separations of these systems indicates 
that if binary formation was dominated by the L3 
mechanism, the velocity dispersion in the primor- 
dial disk must have been less than a few times 
the Hill velocity. A low velocity dispersion sup- 
ports the observed preference for low mutual incli- 
nations. However, the velocity dispersion could not 
have been much below the Hill velocity or the L2S 
mechanism would have dominated and the ultra- 
wide binaries should exhibit a strong preference for 
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retrograde orbits. Thus, under the assumption that 
these are the only two efficient formation pathways, 
it appears that v vh ■ Further modeling of these 
mechanisms is required to determine whether a bal- 
ance can be struck between them where widely- 
separated binaries can be formed with aligned poles 
and roughly equal numbers of prograde and retro- 
grade orientations. 

6. Current simulations of the alternative formation 
mechanism of gravitational collapse and fragmen- 
tation (Nesvorny et al. 2010) create orbital distri- 
butions similar to the observed present-day distri- 
butions. These simulations tend to over-produce 
wide systems, which leaves some room for post- 
formation disruption of ultra-wide binaries. How- 
ever, it remains unclear if this mechanism can pro- 
duce a large number of retrograde systems. 

7. Assuming realistic densities, the implied albedos 
for the ultra-wide TNBs range over ^ 0.09 — 0.30, 
consistent with estimates of the albedos of solitary 
Cold Classical objects (eg., Brucker et al. 2009). 
The distribution of observed albedos with primary 
radius for all dynamically-cold TNBs (seven from 
this work and four others from literature) suggests 
that albedos at the lower end of this scale are in- 
trinsically more common in this population. A 
Gaussian albedo distribution, centered at p = 0.05 
and clipped such that p > 0.05, is consistent with 
observations for widths 0.058 < (Jp < 0.1. This 
estimate does not account for the possibility that 
binary fraction might vary substantially with pri- 
mary radius as predicted by Nesvorny et al. (2011). 

8. All seven systems characterized in this work are 
widely separated enough to have a substantial 
probability of disruption if ever subjected to close 
encounters with Neptune (Parker & Kavelaars 
2010), suggestive of in situ formation or a more 
gentle migration mechanism. Additionally, the col- 
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lisional lifetimes of these binary systems are short 
if the size distribution of the Kuiper Belt is steep 
{q ~ 3.5) at small sizes; for the present population 
of binaries to have survived for the age of the solar 
system requires that they were never subjected to 
a period of intense collisional grinding (eg.. Petit 
& Mousis 2004, Nesvorny et al. 2011) and that 
significant collisions have remained relatively rare 
in the Kuiper Belt over the age of the solar sys- 
tem. We explore the collisional evolution of these 
binary systems and further constrain the collisional 
environment of the present-day Kuiper Belt in an 
upcoming paper. 

9. A confluence of results — (a) gravitational col- 
lapse producing similar binary properties as are 
observed, (b) the evidence that the cold Classical 
Kuiper Belt was not subjected to a period of close 
encounters with Neptune, and (c) the mutual or- 
bits of the wide binaries suggesting formation in 
a very dynamically-cold environment, yet the L2S 
mechanism not dominating the binary formation 
process — all suggest that the cold Classical Kuiper 
Belt formed in situ through rapid accretion of small 
particles directly into large (10 — 100 km) planetes- 
imals, driven by turbulent concentration of solids 
and gravitational collapse. 
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